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Abstract 

This  paper  introduces  a  semi-empirical  dynamic  model  for  PEM  fuel  cell  power  modules  with  potential  applications  in  distributed  generation. 
The  proposed  model  is  constructed  based  on  the  measurements  from  a  Nexa™  PEM  fuel  cell  power  module  under  different  load  conditions,  and 
the  model  has  been  validated  by  static  as  well  as  dynamic  tests.  The  effects  of  temperature  and  variations  in  the  internal  impedance  under  different 
load  conditions  have  been  studied.  The  results  have  indicated  that  the  model  provides  an  accurate  representation  of  the  dynamic  and  static  behaviors 
of  the  fuel  cell  power  module. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

As  a  clean  energy  conversion  device,  fuel  cells  have  recently 
attracted  a  great  deal  of  attention  with  some  promising  results 
in  applications  ranging  from  powering  small  cellular  phones  to 
large  power  generation  in  utilities  [1].  There  are  several  dif¬ 
ferent  types  of  fuel  cells  depending  on  the  type  of  electrolyte 
materials  used  [2] .  Each  fuel  cell  type  has  its  own  characteris¬ 
tics.  To  design  control  systems  for  fuel  cells  and  to  compensate 
its  dynamic  interactions  with  rest  of  the  system,  it  is  necessary 
to  have  a  reasonably  accurate  dynamic  model.  In  fact,  there  are 
several  fuel  cell  models  reported  in  the  literature;  some  are  even 
commercially  available  [3].  Generally  speaking,  the  models  can 
be  categorized  into  two  types,  theoretical  models  based  on  phys¬ 
ical  conservation  laws  [4-12]  and  semi-empirical  models  based 
on  experiments  [13,14].  Different  models  have  been  developed 
for  different  type  of  fuel  cells.  For  example,  the  model  proposed 
in  [5]  is  for  molten  carbonate  fuel  cells  (MCFC)  for  power  gen¬ 
eration;  the  models  discussed  in  [7,9]  are  for  solid  oxide  fuel 
cells  (SOFC).  There  are  also  many  papers  discussing  the  models 
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for  proton  exchange  membrane  (PEM)  fuel  cells  [4,6,8,10-14]. 
PEM  fuel  cells  generally  operate  at  lower  pressure  and  lower 
temperature  with  higher  power  density  compared  to  other  types 
of  fuel  cells.  Therefore,  they  are  more  suitable  for  applications 
in  small  to  medium  power  levels,  such  as  fuel  cell  powered  auto¬ 
mobiles  or  micro-grid  power  applications. 

To  design  an  adequate  fuel  cell  based  distributed  power  gen¬ 
eration  system  to  accommodate  different  load  changes,  it  is 
essential  to  have  an  accurate  dynamic  model  for  the  fuel  cell 
system  so  that  adequate  control  systems  can  be  designed  to  meet 
the  load  demand.  In  such  a  situation,  the  dynamics  of  the  temper¬ 
ature  and  the  internal  resistance  characters  tics  of  fuel  cells  have 
to  be  considered.  Padulles  et  al.  introduces  a  detailed  model  for 
the  SOFC  simulation  in  [7];  and  El-Sharkh  et  al.  subsequently 
modifies  this  model  for  the  simulation  of  PEM  fuel  cells  in  [8]. 
However,  both  have  assumed  the  temperature  and  the  internal 
resistance  to  be  constant.  In  fact,  both  the  stack  temperature  and 
internal  resistance  will  have  effects  on  the  static  and  dynamic 
characteristics  of  the  fuel  cell  output.  Although  such  effects  may 
not  be  so  significant  in  the  case  of  large  scale  fuel  cells,  where  the 
large  thermal  mass  may  prevent  an  excessive  temperature  excur¬ 
sion,  however,  for  smaller  fuel  cells  in  the  range  of  a  few  kW 
levels,  such  as  the  Nexa™  PEM  fuel  cell  power  module,  it  has 
been  observed  that  both  the  stack  temperature  and  the  equivalent 
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internal  resistance  of  the  power  module  can  change  dramatically 
with  different  load  conditions  under  different  hydrogen  pres¬ 
sures,  which  consequently  affect  the  output  characteristics  of 
the  fuel  cell.  This  is  particularly  true  if  the  power  module  is  to 
be  used  with  a  variable  load.  Therefore,  such  effects  cannot  and 
should  not  be  ignored.  This  forms  the  motivation  of  the  work 
reported  herein. 

The  objectives  of  the  paper  are  to  investigate  the  effects  of 
the  temperature  and  the  equivalent  internal  resistance  on  the 
output  characteristics  of  a  PEM  power  module  under  different 
load  conditions,  and  to  develop  a  dynamic  model  for  this  mod¬ 
ule.  The  current  investigation  is  based  on  the  work  reported  in 
[8]  by  taking  into  account  of  the  electrochemical  and  thermal 
aspects  of  chemical  reactions  within  the  stack  as  well  as  the 
equivalent  internal  resistance.  The  developed  model  has  been 
validated  by  experiments  on  a  Nexa™  PEM  fuel  cell  power 
module.  The  results  have  confirmed  that  the  dynamic  model  can 
indeed  represent  the  static  and  dynamic  characteristics  of  the 
fuel  cell  system. 

This  paper  is  organized  as  follows.  In  Section  2,  the  Nexa™ 
PEM  power  module  will  be  briefly  described.  Section  3  provides 
some  experimental  details.  The  dynamic  model  constructions 
have  been  presented  in  Section  4.  The  model  has  been  validated 
in  Section  5  followed  by  the  conclusions. 

2.  Nexa™  PEM  fuel  cell  power  module 

PEM  fuel  cell  power  modules  have  many  unique  features 
as  compared  with  other  fuel  cell  types,  such  as  relatively  low 
operating  temperature,  high  power  density,  and  high  modular¬ 
ity.  They  can  be  tailored  to  different  applications,  in  particu¬ 
lar  for  mobile  applications  and  small-scale  power  generation 
[1].  The  power  generated  by  fuel  cells  can  provide  a  viable 
alternative  in  areas  where  noise,  vibration  or  emissions  are  of 
concerns.  Although  it  is  still  more  expensive  in  comparison 
with  bulk  power  generation,  with  advances  in  technologies  and 
reduction  in  the  cost  of  producing  hydrogen,  PEM  fuel  cells 
can,  one  day,  become  more  competitive  and  popular  in  various 
applications. 

The  Nexa™  PEM  power  module  from  Ballard  Power  Sys¬ 
tems  Inc.  is  capable  of  providing  1.2  kW  of  unregulated  dc 
output.  The  output  voltage  level  can  vary  from  43  V  at  no  load 
to  about  26  V  at  the  full  load.  The  designed  operating  tem¬ 
perature  in  the  stack  is  around  65  °C  at  the  full  load.  There 
are  totally  47  cells  connected  in  series  in  the  stack.  The  fuel 
is  99.99%  hydrogen  with  no  humidification,  and  the  hydrogen 
pressure  to  the  stack  is  normally  maintained  at  5  psig.  Oxygen 
comes  from  the  ambient  air.  The  pressure  of  the  oxidant  air 
is  2.2  psig,  and  the  air  is  humidified  through  a  built-in  humid¬ 
ity  exchanger  to  maintain  membrane  saturation  and  prolong  the 
life  of  the  membrane.  There  is  a  small  compressor  supplying 
excess  oxidant  air  to  the  fuel  cell,  and  the  speed  of  the  com¬ 
pressor  can  be  adjusted  to  match  the  power  demand  from  the 
fuel  cell  stack.  The  Nexa™  fuel  cell  stack  is  air-cooled;  the 
cooling  fan  draws  air  from  the  ambient  surroundings  in  order  to 
cool  the  fuel  cell  stack  and  regulate  the  operating  temperature 
[15]. 


After  start-up,  the  power  module  will  provide  all  the  neces¬ 
sary  internal  power  requirements  for  operation,  such  as  cooling 
fan,  air  compressor,  and  the  control  circuit  board.  Therefore,  the 
output  characteristics  of  the  power  module  can  be  somewhat 
different  from  that  of  the  fuel  cell  stack  itself.  Since  the  power 
module  will  be  used  as  a  unit  and  be  integrated  into  a  distributed 
power  system,  it  is  necessary  to  have  a  complete  knowledge  of 
the  dynamic  characteristics  of  the  entire  power  module,  rather 
than  just  that  of  the  fuel  cell  stack. 

3.  Experimental  investigation  of  the  temperature  profile 
and  the  equivalent  internal  resistance 

3.1.  Experimental  setup 

Through  experiments,  the  characteristics  of  the  stack  tem¬ 
perature  and  the  equivalent  internal  resistance  of  the  Nexa™ 
power  module  have  been  investigated.  To  study  the  effect  of 
the  change  in  hydrogen  pressures  on  the  output  characteristics 
of  the  power  module,  a  hydrogen  bypass  with  an  adjustable 
pressure  regulator  has  been  installed.  The  test  bed  is  shown  in 
Fig.  1. 

3.2.  Stack  temperature  characteristics 

To  investigate  the  change  in  temperatures  for  different  load 
conditions,  the  hydrogen  pressure  is  kept  at  5  psig.  The  load 
connected  to  the  power  module  is  adjusted  by  varying  load  resis¬ 
tances.  The  results  of  these  tests  are  shown  in  Fig.  2.  The  stack 
temperature  reaches  3 1  °C  from  the  room  temperature  (25  °C)  in 
about  5  min  when  5  A  of  the  load  current  is  drawn.  The  tempera¬ 
ture  reaches  44  °C  in  about  5  min  if  the  load  current  is  increased 
to  20  A.  If  the  output  current  is  further  increased  to  47  A,  the  tem¬ 
perature  will  reach  to  66.5  °C  in  about  4  min.  Fig.  2  also  indicates 
that  as  the  load  changes,  the  stack  current  follows  immediately 
with  no  noticeable  time  delays.  Subsequently,  the  temperature 
of  the  fuel  cell  stack  also  changes,  and  the  rate  of  the  change 
is  proportional  to  the  level  of  the  load.  It  should  be  noted  that 
there  are  some  variations  of  the  cell  temperature  along  the  stacks 


Fig.  1.  Experimental  setup. 
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Time  (s) 

Fig.  2.  Temperature  profile  as  a  function  of  stack  current. 


from  the  inlet  to  the  outlet  as  discussed  in  the  literature;  however, 
the  Nexa™  fuel  cell  stack  is  air-cooled,  the  cooling  fan  draws 
air  from  the  ambient  not  from  the  oxidant  stream  so  as  to  cool 
the  fuel  cell  stack  uniformly  [15].  Therefore,  those  variations 
can  be  negligible,  and  the  measured  temperature  is  an  average 
value. 


3.3.  Equivalent  internal  resistance 

As  mentioned  in  Section  2,  after  start-up,  the  Nexa™  power 
module  will  provide  all  necessary  internal  power  requirements 
for  sustaining  the  operation,  such  as  cooling  fan,  air  compressor, 
and  the  control  circuit  board.  These  devices  cannot  be  separated 
from  the  module,  therefore,  they  should  be  considered  as  an 
integral  part  of  the  module.  As  a  result,  the  equivalent  inter¬ 
nal  resistance  R[nt  of  the  module  is  different  from  the  internal 
resistance  of  fuel  cell  stacks  which  is  often  considered  in  the 
literature,  even  though  it  may  be  true  that  the  big  part  of  Rmt 
is  due  to  the  internal  resistance  of  the  stack.  When  the  power 
module  is  used  as  an  energy  source,  it  is  more  accurate  to 
use  R[nt  to  measure  the  internal  resistance.  To  investigate  the 
characteristics  of  this  equivalent  internal  resistance  R{nt,  a  resis¬ 
tive  load  bank  is  used,  and  the  periodic  current  interruption 
(PCI)  mechanism  is  employed  to  maintain  an  isothermal  stack 
temperature  [16].  The  experimental  circuit  is  shown  in  Fig.  3. 
The  0.6  £2  load  resistor  is  frequently  switched  in  and  out  by 
the  switch  S  in  order  to  maintain  the  stack  temperature  at  the 
desired  levels  at  which  the  R[nt  is  measured.  In  Fig.  3,  VQ  is 
the  open-circuit  voltage  of  the  power  module.  7fc  is  the  module 
output  current;  Vl  is  the  output  voltage  of  the  power  module 
when  a  load  is  connected.  The  load  is  adjusted  to  obtain  dif¬ 
ferent  output  voltages  of  the  module.  R[nt  can  be  calculated  by 
(1). 


V^-Vl 

Ifc 


The  output  characteristics  of  the  Nexa™  module  at  differ¬ 
ent  stack  temperatures  and  different  hydrogen  pressures  have 
been  examined.  The  relationship  between  the  equivalent  internal 
resistance  R[n{  and  the  module  output  current  has  been  recorded 
and  illustrated  in  Fig.  4. 

When  a  fuel  cell  is  connected  to  a  load,  the  output  volt¬ 
age  is  always  less  than  the  open-circuit  voltage,  because  there 
are  some  voltage  drops  across  the  fuel  cell  caused  by  the  acti¬ 
vation  loss,  ohmic  resistance  voltage  drop,  and  concentration 
overpotential.  Therefore,  there  are  equivalent  activation  resis¬ 
tance,  equivalent  ohmic  resistance,  and  equivalent  concentra¬ 
tion  resistance  correspondingly.  The  total  internal  resistance  of 
fuel  cell  stack  consists  of  these  resistances,  which  are  current- 
dependent  and  temperature-dependent  [17],  and  it  mainly  deter¬ 
mines  the  equivalent  internal  resistance  of  the  power  module. 
Fig.  4  clearly  shows  that  when  the  stack  current  increases, 
the  equivalent  internal  resistance  dramatically  decreases.  This 
is  also  in  part  because  the  conductivity  of  the  membrane 
is  largely  affected  by  the  water  content  and  the  temperature 
[18]. 

The  results  of  the  experiments  have  indicated  that  both  the 
stack  temperature  and  the  equivalent  internal  resistance  change 
dramatically  with  the  change  in  the  load  currents.  Consequently, 
the  assumption  of  a  constant  stack  temperature  and  internal  resis¬ 
tance  would  not  be  adequate  in  the  case  of  a  Nexa™  PEM  power 


Fig.  4.  The  relationship  between  Rm{  and  the  module  output  current. 
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module.  Such  changes  have  to  be  taken  into  consideration  when 
constructing  a  dynamic  model. 

3.4.  Dynamic  responses  of  the  power  module 

To  investigate  the  dynamic  response  of  the  Nexa™  power 
module,  the  output  voltage  is  measured  when  the  module  is 
subject  to  load  increases  from  0  per  unit  (pu)  to  0.41,  0.73, 
and  0.88  pu,  respectively,  and  also  to  load  decreases  from  0.41, 
0.73,  and  0.88  to  Opu,  respectively,  under  the  nominal  hydro¬ 
gen  pressure  of  5  psig  (1  pu=  1.2  kW,  i.e.  the  full  capability  of 
the  system).  The  results  are  illustrated  in  Figs.  5  and  6  for  both 
cases.  It  can  be  seen  that  the  output  voltage  changes  almost 
instantaneously  when  the  load  changes,  but  it  takes  more  than 
1 .5  min  for  the  output  to  reach  a  new  steady  state.  These  observa¬ 
tions  suggest  that  some  kind  of  independent  energy  sources  such 
as  batteries  may  be  needed  to  supply  the  instantaneous  power 
demand  so  as  to  improve  the  transient  performance  of  the  overall 
system. 


4.  Dynamic  model  construction 


4.1.  Thermal  dynamics 


Using  a  curve  fitting  technique  with  the  experimental  data 
of  Fig.  2,  the  dynamic  characteristics  of  the  temperature  change 
can  be  approximated  by  Eqs.  (2)  and  (3). 

Tm  =  T0  +  (T0  -  Tn  +  7ic/fc)(l  -  e“f/Tr)  (2) 


where  the  time  constant  can  be  expressed  as: 

T\t 

Jt  =  —  (3) 

1  fc 


and  Tm  is  the  measured  temperature  of  the  fuel  cell  stack  (°C), 
7fc  is  the  stack  current  (A),  To,  Trt,  7jc,  and  Tn  are  the  empirical 
parameters. 

For  this  particular  power  module,  the  best  match  occurs  with 
7o  =  28°C,  Tn  =  20,  7jc  =  0.7,  and  7jt  =  4000  at  different  stack 
currents  (5,  20,  and  47  A).  The  results  are  shown  in  Fig.  7. 

By  taking  the  Laplace  transform,  the  transfer  function  of 
the  temperature  change  with  respect  to  the  load  current  can  be 
expressed  in  Eq.  (4). 


T(s)  =  T0  + 


To  Trt  +  T\cIfc 

Tjs  +  1 


4.2.  Characteristics  of  the  equivalent  internal  resistance 


The  load  characteristics  of  the  module  are  investigated  under 
different  stack  temperatures.  Based  on  the  experimental  results 
as  shown  in  Fig.  4,  the  equivalent  internal  resistance  of  the  power 
module  can  best  be  described  by: 

flint  =  Ar  +  Ro  exp  ~  bR  Wife)  (5) 

where  R[nt  is  the  equivalent  internal  resistance,  Ar,  Z?r,  Rq,  tr 
are  all  empirical  parameters. 


Fig.  7.  Comparison  of  the  measured  and  the  computed  stack  temperatures  from 
the  proposed  model. 
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Isolating  the  hydrogen  partial  pressure  and  taking  the  Laplace 
transform,  the  dynamic  characteristics  of  the  fuel  cell  can  be 
written  as: 

PHo  (s)  =  (<4  -  2KTIfc)  (9) 

1  +  Th£  2 


and 


m2 


Kh2RT 


(10) 


where  Kn2  is  the  hydrogen  valve  molar  constant 
(kmol  (atm  s)-1). 

The  partial  pressure  of  the  other  reactants  and  products  can  be 
derived  in  a  similar  way.  Applying  Nemst’s  equation  and  Ohm’s 
law,  the  output  voltage  of  the  fuel  cell  can  be  expressed  in  Eq. 
(11). 


Fig.  8.  Comparison  of  the  computed  equivalent  internal  resistance  and  the  mea¬ 
sured  R[nt. 


To  achieve  reasonably  accurate  representation,  the  following 
values  have  been  selected:  Ar  =  0.82,  Z?r  =  0.13,  7?o  =  0.8,  and 
tr  =  5.  The  synthesized  curves  and  the  measurement  data  are 
shown  in  Fig.  8.  These  curves  indicate  a  good  match  between  the 
computed  values  based  on  the  proposed  models  and  the  exper¬ 
imental  data.  According  to  Eq.  (5),  a  resistance  dynamic  block 
of  the  Nexa™  module  has  been  developed,  and  integrated  into 
the  proposed  dynamic  model. 


4.3.  Dynamic  model  of  the  fuel  cell  module 


For  hydrogen,  the  derivative  of  the  partial  pressure  can  be 
expressed  as  follows  [7-9]: 


d/?H2 
d  t 


where  R  is  the  universal  gas  constant  (1  atm  (kmol  K)-1);  T  is 
the  absolute  temperature  (K);  Van  is  the  volume  of  the  anode 
(1);  Pu2  is  the  hydrogen  partial  pressure  (atm);  ql^2  is  the  hydro¬ 
gen  input  flow  (kmol  s-1);  qrUn  is  the  hydrogen  flow  that  reacts 
inside  the  stack  (kmol  s-1);  and  q ^  is  the  hydrogen  output  flow 

(kmol  s-1).  In  small  size  PEM  fuel  cell  power  modules  such  as 
the  Nexa™  module,  qff2  can  be  set  to  zero,  because  the  channel 
of  the  anode  is  usually  a  dead-end. 

According  to  the  fundamental  electrochemical  relationship, 
the  reacting  hydrogen  flow  can  be  calculated  as: 


M)/fc 

IF 


where  F  is  the  Faraday’s  constant  (C  (kmol)-1);  KY  is  the  mod¬ 
eling  constant  (kmol  (s  A)-1);  and  No  is  the  number  of  the  cells 
connected  in  series  in  the  stack,  which  is  47  for  the  Nexa™ 
module. 


Vfc  =  N0 


RT 

Eo  H - In 

IF 


(  ph2Po52 
\  Pa2o 


FJ—  in  (— 

2a  F  \  io 


^intAc 


(ID 


Based  on  the  above  equations,  a  dynamic  model  of  the 
Nexa™  power  module,  which  includes  the  effects  of  tempera¬ 
ture  and  internal  resistance  variations,  has  been  developed  and 
validated  with  the  experiments. 


5.  Model  validation  and  discussion 


5.1.  Static  output  characteristics 

The  comparison  of  fuel  cell  characteristics  based  on  the 
work  in  [8],  static  output  characteristics  based  on  Eq.  (11), 
and  the  measured  data  have  been  carried  out.  The  results  are 
presented  in  Fig.  9.  The  dashed  line  marked  as  ‘Result  from 
[8]’  is  from  [8],  which  assumes  a  constant  temperature  and  a 
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(a)  Time  (s) 


Fig.  10.  Comparison  of  dynamic  responses  subject  to  step  load  changes:  (a)  a 
step  load  increase  and  (b)  a  step  load  decrease. 


constant  internal  resistance.  The  solid  line  is  traced  from  the 
experimental  data,  while  the  dotted  line  with  circles  is  com¬ 
puted  from  the  model,  where  effects  of  temperature  and  varia¬ 
tions  in  the  equivalent  internal  resistance  have  been  considered. 
The  results  clearly  show  that  the  developed  model  matches 
significantly  better  than  that  from  the  model  proposed  in  [8] 
with  the  assumptions  of  constant  temperature  and  the  internal 
resistance. 

5.2.  Dynamic  output  characteristics 

To  validate  the  dynamic  characteristics  of  the  output  volt¬ 
age  with  respect  to  load  changes,  a  step  load  increase  from  0  to 
0.73  pu  and  a  step  load  decrease  from  0.73  to  0  pu  have  been  car¬ 
ried  out.  The  results  of  both  the  model  of  [8]  and  the  developed 
model  are  shown  in  Fig.  10.  The  comparisons  of  experimen¬ 
tal  data  with  these  models  are  also  included  in  these  figures. 
Clearly,  a  good  agreement  between  the  predicted  response  from 


the  developed  model  and  the  measured  data  has  been  established. 
The  results  also  indicate  that  there  is  a  significant  discrep¬ 
ancy  between  the  model  proposed  in  [8]  and  the  experimental 
data. 

6.  Conclusions 

Through  experiments,  this  paper  conclusively  show  that  both 
the  stack  temperature  and  the  equivalent  internal  resistance  of  the 
Nexa™  PEM  power  module  have  significant  effects  on  the  static 
as  well  as  dynamic  characteristics  of  the  fuel  cell  output.  Conse¬ 
quently,  they  should  not  be  ignored.  A  semi-empirical  dynamic 
model  for  a  Nexa™  PEM  power  module  considering  the  charac¬ 
teristics  of  the  temperature  and  the  equivalent  internal  resistance 
has  been  developed  in  the  paper.  This  developed  model  has  been 
validated  using  the  experimental  data.  The  results  have  convinc¬ 
ingly  shown  that  the  developed  model  can  accurately  represent 
the  Nexa™  PEM  power  module  in  a  wide  range  of  load  condi¬ 
tions.  Consequently,  this  model  is  more  suitable  for  applications 
in  which  the  load  is  not  a  constant. 
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